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Abstract We report here on X-ray and IR observations 
of the Anomalous X-ray Pulsar (AXP) IRXS J170849- 
400910. First, we report on new XMM-Newton, Swift- 
XRT and Chandra observations of this AXP, which con- 
firm the intensity-hardness correlation observed in the 
long term X-ray monitoring of this source. These new X- 
ray observations show that the AXP flux is rising again, 
and the spectrum hardening. If the increase of the source 
intensity is indeed connected with the glitches and a pos- 
sible bursting activity, we expect this source to enter in a 
bursting active phase around 2006-2007. Second, we re- 
port on deep IR observations of IRXS J170849-400910, 
taken with the VLT-NACO adaptive optics, showing that 
there are many weak sources consistent with the AXP 
position. Neither star A or B, as previously proposed by 
different authors, might yet be conclusively recognised 
as the IR counterpart of IRXS J170849-400910. Third, 
using Monte Carlo simulations, we re-address the calcu- 
lation of the significance of the absorption line found in a 
phase-resolved spectrum of this source, and interpreted 
as a resonant scattering cyclotron feature. 
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1 Introduction 

AXPs are a small group of neutron stars (NSs) which 
stand apart from other known classes of X-ray sources. 
At the moment there are 7 confirmed AXPs plus 2 can- 
didates. These X-ray pulsars share many peculiarities: 
they are all (but one) radio-quiet (Camilo et al. 2006; 
Burgay et al. this volume), exhibit X-ray pulsations with 
spin periods in a small range of values (~5-12s), they 
have a large spin-down rate {P « 10~^^ — 10~^°ss~-^), a 
rather high X-ray luminosity {Lx ~ 10^'' — lO^^ergs""'^), 
and faint IR counterparts with i^s ~ 20 — 22 magnitudes 
(for a recent review see Woods & Thompson 2004 and 
Kaspi 2006 in this volume). The nature of their X-ray 
emission was intriguing all along. In fact, it is too high 
to be produced by the loss of rotational energy alone, but 
on the other hand, no hints for a companion star were 
found, neither through deep observations at other wave- 
length, nor timing the X-ray pulsations with the hope 
of finding Doppler shifts (Israel et al. 2003a; Mereghetti, 
Israel & Stella 1998). 

At present, the model which is most successful in ex- 
plaining the peculiar observational properties of AXPs is 
the "magnetar" model. In this scenario AXPs are thought 
to be isolated NSs endowed with ultra-high magnetic 
fields {B ~ lO^"* — 10^^ Gauss) and their X-ray emis- 
sion powered by magnetic field decay (Duncan & Thomp- 
son 1992; Thompson & Duncan 1993, 1996). This idea 
is strongly supported by the estimate of AXPs' mag- 
netic field through the classical dipole braking formula, 

B ~ 3.2 X 10^^ V Pp Gauss, which gives in all cases val- 
ues above the electron critical magnetic field {Bqed ~ 
4.4 X 10^^ Gauss). 
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Table 1 Best fit values of the spectral parameters obtained for about ten years X-ray monitoring of IRXS J170849-400910 . 
Fluxes (and percentages of fluxes) are unabsorbed, in units of lO"'^" ergs~^ cm~^ and in the 0.5-10 keV energy range. The 



Nh was fixed at the XMM-Newton value of 1.36 x 10^ 
confidence level. 



for all the observations. See text for details. Errors are at 90% 



ASCA1996 SAX 1999 SAX 2001 Chandra 2002 XMM2003 Chandra 2004 Swift 2005 



kT (keV) 

Gamma 

Flux 

PL Flux (%) 
Xg (d.o.f.) 



41+0 01 
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+0.03 
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83 ±3 



2.77 
1.30 



0.430^ 



"1+0.015 
-0.017 

2.62+°°* 



0.02 
+0.01 

0.04 

71 ±3 



1.43 



0.95 (147) 1.11 (182) 



Alternative scenarios, invoking accretion from a fossil 
disk remnant of the supernova explosion (van Paradijs, 
Taam & van den Heuvel 1995; Chatterjee, Hernquist & 
Narayan 2000; Alpar 2001), are still open possibilities 
although encounter increasing difficulties in explaining 
the data. 

1RXSJ170849-400910 was discovered with ROSAT 
(Voges et al. 1996) and later on a ^^lls modulation 
was found in its X-ray flux with ASCA (Sugizaki et al. 
1997). Early measurements suggested that it was a fairly 
stable rotator with a spin period derivative of ~ 1.9 x 
10~^^ss~^ (Israel et al. 1999). However, in the last four 
years the source experienced two glitches, with differ- 
ent post-glitch recoveries (Kaspi, Lackey & Chakrabarty 
2000, DairOsso et al. 2003, Kaspi & Gavriil 2003). Searches 
for optical/IR counterparts ruled out the presence of a 
massive companion (Israel et al. 1999). Very recently, two 
different objects were proposed by different groups, as 
being IRXS J170849-400910 IR counterpart, and there 
still is an open debate on which one is the AXP counter- 
part (Israel et al. 2003b; Safi-Harb & West 2005; Durant 
& van Kerkwijk 2006). A diffuse (^ 8') radio emission 
at 1.4 GHz was recently reported, possibly associated 
with the supernova remnant G346. 5-0.1 (Gaensler et al. 
2000). 

Pulse phase spectroscopy analysis of two i?eppoSAX 
observations of IRXS J170849-400910 (Israel et al. 2001; 
Rea et al. 2003) revealed i) a large spectral variability 
with the spin-phase, ii) a strong energy dependence of 
the pulse profile shape, and iii) shifts in the pulse phase 
between the low and the high energy profiles. High vari- 
ability of the pulse shape with energy is now detected at 
higher energies, up to ~220 keV (Kuiper et al. 2006). 

By analysing a iJeppoSAX observation taken in 2001 
(the longest pointing ever performed on this source), Rea 
et al. (2003) reported the presence of an absorption line 
at ~ 8 keV in a phase-resolved spectrum. Interpreting 
the feature as a cyclotron line due to resonant scattering 
yields a neutron star magnetic field of either 9.2 x 10^"'^ G 
or 1.6 X 10^^ G, in the case of electron or proton scatter- 
ing, respectively. 

In Sllwe report on new IRXS J170849-400910 XMM- 
Newton, Chandra and Swift observations which were used 
together with the previous ROSAT, ASCA and BeppoSAX 
observations to monitor the X-ray spectrum and flux of 



the AXP. Then we report on deep infrared VLT-NACO 
observations of this AXP. In ijSlwe carefully re-address 
the calculation of the significance of the absorption line 
found around 8keV during a long BeppoSAX. observa- 
tion (Rea et al. 2003), we then summarise and discuss 
all the results in fQ] 



2 Observations and Results 

In this section we report on the observations and the data 
analysis of the four new X-ray and the IR observations. 
The results on the X-ray timing analysis are reported 
below in the text while the X-ray spectral parameters 
may be found in Tab. 1. All the X-ray spectra were fit 
by an absorbed blackbody plus a power-law component 
(but see also Rea, Zane, Lyutikov & TuroUa in this vol- 
ume for a different spectral modelling). Giving the very 
high statistics we have in the XMM-Newton observa- 
tion, and assuming the absorption does not vary along 
the line of sight, we fixed for all the spectra the absorp- 
tion at the XMM-Newton value (phabs XSPEC model: 
Nh = (1.36±0.04) x 10^^ cm~^; abundances from Anders 
& Grevesse 1989) 

2.1 XMM-Newton 

IRXS J170849-400910 was observed with XMM-Newton 
between 2003 August 28th and 29th, for - 50 ks. The 
MOS cameras were operated in Prime Partial Window 
Mode, while the PN camera was in Prime Small Window 
Mode, all with the medium optical photons blocking fil- 
ter. Since a higher background affected the last '^lOks of 
the observation, we used only the data during intervals 
in which the count rate above 10 keV was less than 0.35 
counts s~^. The source events and spectra were extracted 
within a circular region of 27" centred on the peak of the 
point spread function of the source. This non standard 
radius was used because the source was located near the 
edge of the chip. The background was obtained from a 
source-free region of 27". In order to determine the spin 
period of IRXS J170849-400910 we barycentered the 
events arrival times and obtained, through a phasefitting 
technique, a best spin period of Ps = 11. 00170±0. 00004s 
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Fig. 1 VLT-NACO image in the Ks band of the field of 
IRXS J170849-400910 . The Chandra 0.8" 90% error circle 
is over-plotted (Israel et al. 2003b). We marked all the faint 
sources being consistent at 3(j with the AXP position. Note 
that our astrometry (following Israel et al. 2003b) is slightly 
different from the one reported by Durant & van Kerk- 
wijk (2006), this is due to a different catalogue used for the 
astrometry. 

(all errors are at the 90% confidence level). We found 
that the pulsed fraction of the X-ray signal (defined as 
the amplitude of the best-fitting sine wave divided by 
the, background corrected, constant level of the emis- 
sion) is energy-dependent, and it varies from 39.0 ±0.5% 
in the 0.5-2.0 keV range to 29 ±1.5% in the 6.0-10.0 keV 
range. These values are consistent with those reported 
for the pre-glitches iJeppoSAX observation (Israel et al. 
2001) while both are larger than those reported for the 
post-glitches BeppoSAX observation (Rea et al. 2003). 
Detailed results for this observation are reported in Rea 
et al. (2005a). 

2.2 Chandra 

IRXS J170849-400910 was observed by the Chandra M- 
vanced CCD Imaging Spectrometer (ACIS), first for ~30 ks 
with the High Energy Transmission Grating Spectrome- 
ter (HETGS) on 2002 September 9th, then for ~30ks in 
Continuous Clocking (CC) mode on 2004 July 3. For a 
more detailed description of the instruments and on the 
data processing we defer to the Chandra X-ray Center 
(CXC) documents^. Detailed results for these two obser- 
vations are reported in Rea et al. (2005a) and Campana 
et al. (2006). 

2.2.1 High Energy Transmission Crating Spectrometer 



Medium Energy Gratings (MEGs) with a range of 2.5- 
31 A (0.4-5.0 keV) and the High Energy Gratings (HEGs) 
with a range of 1.2-15 A (0.8-10.0 keV). The HETGS 
spectra were imaged by ACIS-S, an array of 6 CCD de- 
tectors normally read-out every 3.2 s. The HETGS/ACIS- 
S combination provides an undispersed (zeroth order) 
image and dispersed spectra from the gratings. The var- 
ious orders overlap and are sorted using the intrinsic en- 
ergy resolution of the ACIS CCDs: A\ = 0.012 A for the 
HEG and 0.023 A for the MEG. 

The MEG and HEG first order count rate were only 
0.5 and 0.2 cts s~-^, we therefore did not expect the 
dispersed spectrum to be affected by pileup, while the 
zeroth-order image was not used in our spectral anal- 
ysis because highly affected by photon pileup. We used 
the standard CIAO tools to create detector response files 
for the MEG and HEG +1 and —1 order spectra. These 
were combined when the +/— order spectra were added 
for the HEG and MEG separately. We binned the data at 
0.08 A with a minimum of 30 counts per bin. To look for 
high-resolution spectral features, the data were binned at 
0.015 A for the HEG and 0.03 A for the MEG. We also 
created background files for the HEG and MEG spectra 
using the standard CIAO tools. 

2.2.2 Continuous Clocking 

In order to avoid pile-up effects, in this second observa- 
tion the source was observed in the Continuous Clock- 
ing (CC) mode (CC33_FAINT; time resolution 2.85ms). 
The source was positioned in the back-illuminated ACIS- 
S3 CCD on the nominal target position. A detailed de- 
scription on the analysis procedures, such as extraction 
regions, corrections and filtering applied to the source 
events and spectra can be found in Rea et al. (2005b). 

In order to perform a timing analysis we corrected 
the events arrival times for the barycenter of the so- 
lar system (with the CIAO axbary tool) using the pro- 
vided ephemeris. We used for the timing analysis only 
the events in the 0.3-8 keV energy range and the stan- 
dard Xronos tools (version 5.19). One fundamental peak 
plus one harmonic were present in the power-spectrum. 
A period of 11.00223 ± 0.00005 s was detected referred 
to MJD 53189. The pulse profile did not change with re- 
spect to the previous detection and the 0.3-8 keV PF is 
35.4 ± 0.5%. Being the CC mode not yet spectrally cal- 
ibrated, the Timed Exposure (TE) mode response ma- 
trices (rmf) and ancillary files (arf) are generally used 
for the spectral analysis. We defer to Rea et al. (2005b) 
for a detailed description on the extraction procedures 
of the spectral matrices. 



The High Energy Transmission Grating Spectrometer 
(HETGS) employs two sets of transmission gratings: the 

^ http://asc.harvard.edu/ud ocs/docs/docs.html: 
|http://asc.harvard.edu/ciao/| 



2.3 Swift 

IRXS J170849-400910was observed with the Swift satel- 
lite a few times between 2005 January 29th and March 
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29th, being a calibration source for the timing accuracy 
and for the wings of the Point Spread Function of the 
X-Ray Telescope (XRT). Here we focus on data taken 
in Window Timing (WT) and Photon Counting (PC) 
mode longer than 1 ks. We extracted data from two 
WT observations. The extraction region is computed au- 
tomatically by the analysis software and is a box 40 
pixels along the WT strip, centred on source, encom- 
passing ~ 98% of the Point Spread Function in this 
observing mode. We extracted photons from PC data 
from an annular region (3 pixels inner radius, 30 pixel 
outer radius) in order to avoid pile-up contamination. 
We consider standard grades 0-2 in WT and 0-12 in PC 
modes. Background spectra were taken from close-by re- 
gions free of sources. The photon arrival times were cor- 
rected to the Solar system barycenter. A period search 
led to a clear detection of the neutron star spin period at 
P = 11. 0027 ± 0.0003 s, derived with a phase fitting tech- 
niques. This period is consistent with the extrapolation 
from known ephemerides at a constant period derivative 
(Kaspi & Gavriil 2003; Dall'Osso et al. 2003). We found 
a PF of 31±2%, 39±3%, 29±4% and 35±7% in the 0.2- 
10 keV, 0.2-2 keV, 2-4 keV and 4-10 keV energy bands, 
respectively. Detailed results are reported in Campana 
et al. (2006). 



2.4 IR observation: VLT-NACO 

A deep observation of the IRXS J170849-400910 field, 
was taken on 2003 June 20th from the Very Large Tele- 
scope using the NAOS-CONICA adaptive optics. We 
defer to Israel et al. 2004 for details in the data reduc- 
tion. In Fig^we present the Kg band field around the 
Chandra 0.8" position of IRXS J170849-400910 (Israel 
et al. 2004). Besides sources A and B proposed by Is- 
rael et al. (2003b), Safi-Harb & West (2005) and Durant 
& van Kerkwijk (2006), as the possible IR counterparts 
to this AXP, many further faint sources (e.g. Star C and 
Star D) were detected strongly consistent with the Chan- 
dra uncertainty region. Unfortunately, H band images 
were obtained under poor sky conditions, and thus the 
H — Ks colour could not be determined for these other 
faint objects. Nonetheless, we note that the Ks magni- 
tudes of e.g. stars C and D (20.3 ± 0.2 and 21.7 ± 0.3, 
respectively) are in better agreement with the IR mag- 
nitudes typical of AXPs, than star A and B previously 
proposed (17.61 ± 0.07 and 18.78 ± 0.05, respectively). 
Deep images in the L' band were also obtained, but no 
object was detected within the Chandra uncertainty re- 
gion at a limiting magnitude of L' ^ 17.8, the deepest 
limit ever obtained for an AXPs in this band. Detailed re- 
sults and discussion will be reported in Israel et al. (2006, 
in prep.). 
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Fig. 2 Unfolded phase-resolved spectrum presenting the ab- 
sorption line discovered in IRXS J170849-400910 (Rea et al. 
2003; Rea et al. 2004). 



3 On the absorption line at 8.1 keV 

During the post glitches BeppoSAX observation in 2001, 
evidence for an absorption feature was found (see Fig|2Jl. 
This feature was not detected during the XMM obser- 
vation 3 years later, the only observation by now with 
a comparable statistics. The XMM-Newton upper limit 
for the line depth is 0.15 at 95% confidence level, which 
compared with the value found by Rea et al. (2003; 
0.8±0.4 at 90% confidence level), leave only a very small 
chance that the two measurements are consistent. 

We then undertook a careful re-analysis of the BeppoSAX. 
data. This re-analysis resulted in the finding that the 
phases at which the absorption line was strongest were 
given incorrectly in the published version of Rea et al. 
(2003), as we noticed earlier: in particular, the line is 
strongest close to the pulse minimum in the 0.1-2 keV 
band (or the pulse maximum in the 6-10 keV band). We 
nevertheless found that the reported estimate of the sig- 
nificance is sound and not much influenced by different 
choices in the background subtraction (annular regions, 
circular regions far from the source or using blank field 
files) or by different extraction regions for the source. 
The re-analysis of the BeppoSAX data made varying 
the extraction radius, the criterion for the background 
subtraction and the spectral binning factor, results in 
basically the same line properties, which strengthen our 
confidence in the robustness of the result. 

Using an F-test method and taking into account the 
six trials we made in the phase resolved spectra, we de- 
rive a confidence level for the absorption line of ~ 4(t. 
Note that even if we take into account all the possible 
energies at which the feature could lie in the LECS plus 
MECS energy range, the confidence level is still <L 3.5cr. 
However, despite being the most common method in 
astrophysics to derive the significance of the emission 
or absorption spectral features, Protassov et al. (2002) 
pointed out that the F-test may be inappropriate in 
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these circumstances, leading sometimes to incorrect sig- 
nificance estimates. 

Following the recipe of Protassov et al. (2002), in or- 
der to further investigate on this significance issue, we 
ran a Monte Carlo simulation of 10"* spectra fixing only 
the continuum model (parameters reported in Rea et al. 
2003) and the same number of photons of the phase 
resolved spectrum which showed the line in the 2001 
BeppoSAK observation. The results of the simulation 
is shown in Fig|3 in the left panel each red circle repre- 
sents one of the 10^ simulated spectra for which a line 
was detected, here we plot the depth of the lines as a 
function of the line energy. On the other hand, in the 
right panel we report the the number of spectra, among 
the 10* simulated spectra, for which the statistical fluc- 
tuation reported the presence of a line at 8.1 keV, as 
the function of the line depth. From this simulation we 
found 32 spectra with depth >0.8 in lO"* points. We can 
then reliably say that the probability of the line being a 
fluctuation is <0.32%. In summary, we confirm the de- 
tection of the 8.1 keV absorption line in the BeppoSAX 
data made by Rea et al. (2003) at 99.68% confidence 
level (see Figl^j). Note that the non homogeneous cover- 
age of red circles over the entire spectral range mirrors 
the energy dependency of the BeppoSAX spectral ma- 
trices. 

The interpretation of the absorption feature as a cy- 
clotron scattering line proposed by Rea et al. (2003) 
was based on the following criteria: 1) a Gaussian line 
gives a bad fit and does not reproduce the asymmetrical 
shape of the observed feature; 2) the best fitting model is 
the XSPEC cyclabs model; 3) the line strength is highly 
phase dependent; 4) no atomic edges or absorption lines 
are known to lie around 8.1 keV (at least without assum- 
ing ad hoc shifts possibly due to the high gravitational 
redshift or to Zeeman effects in such strong magnetic 



field); 5) the relation between the line energy and width 
agrees with that of cyclotron scattering features discov- 
ered in other classes of sources (see Fig. 5 in Rea et al. 
2003); 6) the magnetic field inferred from the line energy, 
either in the case of an electron or proton cyclotron res- 
onance, is reasonably consistent with what is expected 
for a normal neutron star (~ 10^^ G) or for a magnetar 

10^^ G), both being still open possibilities. Then, if 
this feature is real, all the above points hint toward the 
cyclotron nature of the absorption feature at 8.1kev. 

Keeping always in mind the possibility that the ab- 
sorption line in the BeppoSAX spectrum might be due 
to statistical fluctuations, in Rea et al. (2005a) we dis- 
cuss physical mechanisms which could be responsible for 
the appearance of a transient cyclotron line in this source 
in the context of the magnetar scenario. 



4 Discussion 

1. Intensity Hardness correlation 

The long-term evolution of the source flux and the 
spectral hardness are shown in Fig^ where different 
observations, spanning nearly ten years, were collected. 
Comparing the two panels, there seems to be a correla- 
tion between the photon index the source X-ray flux. The 
spectrum became progressively harder as the flux rose 
in correspondence of the two glitches and then softened 
as the luminosity dropped, following the glitch recovery. 
This is suggestive of a scenario in which the mechanism 
responsible for the glitches is also at the basis of the 
enhanced emission and of the spectral hardening. 

The similarity of the second glitch of IRXS J170849- 
400910 with the one discovered during the bursting ac- 
tivity of IE 2259-1-586 (Kaspi et al. 2003), after which 
a similar exponential recovery was seen, suggests that 
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Fig. 4 Correlated photon index and X-ray fiux clianges witli time. Dasiied line represents the two glitch epochs. All fluxes 
are unabsorbed and calculated in the 0.5-10 keV energy band keeping fixed the Nh ~ 1.36 x lO^'^ cm~^. 



bursts likely occurred in IRXS J170849-400910 as well, 
but the sparse observations did miss them, as already 
suggested by Kaspi & Gavriil (2003). Moreover, the spec- 
tral parameters and the flux changes after the recov- 
ery of the glitch strengthens this idea in comparison 
with what was reported for the post-bursts fading of 
IE 2259-1-586 (Woods et al. 2004). 

This may be interpreted as the onset of a twist, which 
grew, culminated in the glitches, and then decayed. A 
twisted external field, in fact, is in an unstable magne- 
tostatic equilibrium and evolves towards a pure dipole 
field which represents the configuration of minimal en- 
ergy (see Rea et al. 2005a for a detailed interpretation). 

The last Chandra and Swift observations show that 
the source is slowly increasing again its flux and harden- 
ing its emission. Whether the suggested correlation with 
the glitching and possibly bursting activity holds, we ex- 
pect the source to re-enter an active state around 2006- 
2007. If confirmed for all the magnetars, the X-ray mon- 
itoring might be an excellent tool to foresee the activity 
of magnetars. 

2. On 1RXSJ170849 400910 IR counterpart 

Based on the deep VLT-NACO observation, we be- 
lieve that the identification of the IR counterpart of IRXS 
J170849 -400910 is still an open issue, mainly due to the 
very crowded region in which this source is located. Note 
that the 2.5ct variability recently reported for source B 
by Durant & van Kerkwijk (2006) does not seem to be 
a conclusive word on the IR counterpart of this source, 
especially considering that source B has IR magnitudes 
much brighter than all other AXPs. 

3. On the possible cyclotron line 

Following Protassov et al. (2002) we re-addressed the 
issue of the significance of the absorption line discovered 



by BeppoSAX . We end up with a 99.68% confidence level 
for its existence. The fact that the source was not com- 
pletely recovered by the second glitch at the time the line 
appeared, make the correlation between line appearance, 
glitching activity and fiux enhancement possibly intrigu- 
ing and might suggests that the conditions for line forma- 
tion were met at the epoch of the BeppoSAX. pointing. 
Future long X-ray observations are needed in order to 
follow the source in its slow flux increase and possibly 
re-detect the source in such state. 
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